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ABSTRACT 

The origin of the so-called '21' /im feature which is especially prominent in 
the spectra of some carbon-rich protoplanetary nebulae (PPNe) is the matter of 
a lively debate. A large number of potential band carriers have been presented 
and discarded within the past decade. The present paper gives an overview of 
the problems related to the hitherto proposed feature identifications, including 
the recently suggested candidate carrier silicon carbide. We also discuss the case 
for spectroscopically promising oxides. 

SiC is shown to produce a strong resonance band at 20-21 yum if coated by 
a layer of silicon dioxide. At low temperatures, core-mantle particles composed 
of SiC and amorphous Si02 indeed have their strongest spectral signature at a 
position of 20.1 /im, which coincides with the position of the '21' ^um emission 
band. 

The optical constants of another candidate carrier that has been relatively 
neglected so far - iron monoxide - are proven to permit a fairly accurate repro- 
duction of the '21' /im feature profile as well, especially when low-temperature 
measurements of the infrared properties of FeO are taken into account. As candi- 
date carrier of the '21' /zm emission band, FeO has the advantage of being stable 
against further oxidation and reduction only in a narrow range of chemical and 
physical conditions, coinciding with the fact that the feature, too, is detected in 
a small group of objects only. However, it is unclear how FeO should form or 
survive particularly in carbon-rich PPNe. 

Subject headings: circumstellar matter — infrared: stars — methods: laboratory — 
stars: AGB and post-AGB — stars: atmospheres 



1. Introduction 

Kwok, Volk & Hrivnak (1989) reported 
the discovery of a strong emission band, 
centered (as it then seemed) at 21 fim, in 
the spectra of four carbon-rich protoplane- 
tary nebuale (PPNe) taken by the Infrared 
Astronomical Satellite (IRAS). The fea- 
ture was not resolved well with the Low 
Resolution Spectrometer (LRS), but nev- 
ertheless, it was evident aheady from the 
LRS observations that it has a consider- 
able strength in some objects (especially 
IRAS 07134+1005) and is too broad (full 
width at half maximum >2 fim) to origi- 
nate from any atomic transition. The pi- 
oneering work by Kwok, Volk & Hrivnak 
also highlighted the fact that the '21' /xm 
feature only occurs in a short transitional 
phase of stellar evolution, i.e. in a very lim- 
ited range of physical and chemical con- 
ditions. Any identification of the '21' /xm 
feature has to cope with this fact. 

Observations of selected PPNe with the 
Infrared Space Observatory (ISO) by Hriv- 
nak, Volk & Kwok (2000) brought sub- 
stantial progress, especially concerning the 
position and width of the feature, which 
was shown to be 20.1 /xm instead of 21 /im 
(notwithstandig, the feature is still being 
refered to as '21' //,m band, a tradition to 
which we also stick in this paper). Hriv- 
nak, Volk & Kwok (2000) analysed ISO 
observations at 2-45 //m for seven PPNe 
and two other carbon-rich objects. Six 
of the sources exhibit emission features at 
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21 and 30 //m. The authors resolved the 

30 /im band into a broad feature peaking 
at 27.2 yum (the '30' fim feature) and a nar- 
rower feature at 25.5 //m (the '26' //m fea- 
ture). The observations suggest that the 
sources with '21' //m feature do also show 
the 30 fim feature. The '26' /xm feature ap- 
pears to be present in a broader range of 
objects - carbon stars, PPNe, and possibly 
PNe - than does the '21' //m feature. 

The candidate carriers of the '21' //m 
band have since its discovery been the mat- 
ter of a long discussion. Cox (1990) pre- 
sented IRAS-LRS-spectra of ten HIl re- 
gions which showed a peak close to 21 //m. 
He associated this peak with the feature 
discovered by Kwok, Volk & Hrivnak and 
tentatively assigned it to iron oxides with 
the stoichiometrics 7-FC2O3 and FC3O4. In 
the late 1990s, some of these H 11 regions 
have also been observed with ISO; in the 
ISO spectra, however, we do not find any 
traces of a broad 21 or 20 ^m emission 
band (in accordance with Oudmaijcr & de 
Winter (1995), who already cast doubt on 
Cox' observations). 

Goebel (1993), still on the basis of IRAS 

spectra, assigned the '21' ^m feature to 
silicon disulfide (SiS2), which has indeed 
an intense vibrational band in the 20 yum 
spectral range (but, unfortunately, a sec- 
ond one which is not observed; see below. 
Sect. 4). Prom a chemical point of view, 
other sulfur compounds like Ss, OCS, S2 
and CS2, even though mainly as mantles 
of other grains, were also taken into ac- 
count by Goebel as potential emitters in 
the 20 //m region. 

Webster (1995) brought fuUeranes as 
potential carriers of the '21' /im feature 
into the discussion. Hill, Jones & d'Hende- 
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court (1998), having studied the IR bands 

of terrestrial nitrogen- doped and neutron- 
irradiated diamonds, showed that defects 
in diamond structures can cause an emis- 
sion band at 20—22 /xm and even stronger 
bands in the 8-10 /xm region. However, 
there is a significant discrepancy between 
the position of 21-22 /im and the actual po- 
sition of the '21' nm band as determined by 
recent observations. 

Papoular (2000) calculated synthetic 
spectra of carbonaceous macromolecules 
containing OH, oxygen, sulfur, nitrogen 
or a combination of them and found that 
it is possible to reproduce not only the 
'21' //m feature, but also the 3.3, 11.3, 
12-I-, 26 and 30 fim bands by emissions 
of such molecules. According to this sce- 
nario, the '21' /im feature would arise from 
out-of-plane vibrations of oxygen or nitro- 
gen atoms in 5-membered carbon rings. 
Papoular's approach seems promising be- 
cause it aims at an identification of all 
the major features detected in the sources 
of the '21' /xm band on the basis of one 
unified dust model. But at the time be- 
ing, substantial problems remain, such as 
the unsatisfactory reproduction of the ob- 
served 26-35 nm band profile. 

There are various other papers dis- 
cussing carbonaceous materials as carriers 
of the '21'//m feature, but without any 
detailed fits of its average profile. Souris- 
seau, Coddens & Papoular (1992) sug- 
gested amides, especially urea, as potential 
carrier. Buss et al. (1990) and Justtanont 
et al. (1996) suggested polycychc aromatic 
hydrocarbon molecules and Grishko et al. 
(2001) hydrogenated amorphous carbon. 

Von Helden et al. (2000) came up 
with another suggestion, attributing the 



20.1 /im feature to titanium carbide clus- 
ters (made of 27 to 125 atoms). Spectra 
of bulk titanium carbide and TiC particles 
as those subsequently published by Hen- 
ning & Mutschke (2001) and Kimura & 
Kaito (2003) show no band in the respec- 
tive wavelength region, but rather weak 
bands at 9.5, 12.5 (Kimura & Kaito 2003) 
and 19 //m (Henning & Mutschke 2001; see 
also Speck & Hofmeister 2004). This does 
not contradict the measurements by von 
Helden et al. (2000), however, since clus- 
ters are likely to have electronic and vi- 
brational properties deviating from those 
of bulk material even if composition and 
crystal structure are similar. 

Notwithstanding, in 2003, three papers 
doubting the 'TiC-cluster-identification' 
with substantial arguments came out, all 
based essentially on arguments concern- 
ing the abundance of titanium (Hony et 
al. 2003, Chigai et al. 2003, and Li 2003). 
Hony et al. (2003) argue that TiC, like any 
kind of dust, can only emit as much in the 
infrared as it absorbs in the UV and visual 
range. Assigning a simple model opacity 
to TiC at short wavelengths, Hony et al. 
found that this energy balance can be ful- 
filled only if unrealistically large amounts 
of TiC are predicted to form in carbon- 
rich PPNe. This point was strengthened 
by Li (2003) by applying a sum-rule based 
on Kramers-Kronig considerations on the 
wavelength-integrated absorptivity. Chi- 
gai et al. (2003) came to the same result, 
comparing the number densities of SiC 
and TiC grains necessary for producing 
the 11 fxm and 21 /xm features, respectively, 
with the available Si and Ti abundances. 

Recently, Speck & Hofmeister (2004) re- 
ported the experimental finding that sili- 
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con carbide, under certain circumstances, 
does not only show the well-known reso- 
nance feature at 11 /im, but also a sec- 
ondary band (not predicted by group the- 
ory), which is centered at 20-21 yum. This 
secondary band is reported to appear only 
in the /9-SiC-polytype, and nitrogen or 
carbon impurities were suspected by the 
authors to favour its occurence. Cold, 
maybe N- or C-doped /3-SiC-dust (with 
temperatures around 100 K) would then, 
according to Speck & Hofmeister, emit 
at 20-21 /im (depending on particle size 
and shape) even more effectively than at 
11 /xm, since the latter peak would still be 
in the Wien regime of the corresponding 
Planck function. Important as this sugges- 
tion might be for the identification of the 
'21' /xm band, it is not beyond any doubt, 
as we try to show in §3. 

The present paper is structured as fol- 
lows. In §2, we summarize the currently 
available astronomical information on the 
sources in which the enigmatic '21' /xm 
band is observed. §3 presents new in- 
sights to the infrared properties of materi- 
als which should, despite some problems, 
not be excluded as feature carriers - such 
as core-mantle particles composed of SiC 
and Si02 and cold FeO dust grains. In §4, 
we deliver additional arguments for ruling 
out potential carriers of this band which 
have been ruled out more or less convinc- 
ingly by other authors. The last section 
(§5) summarizes our conclusions. 

2. Properties of the sources of the 
'21' //m feature 

So far, the '21' //m feature is observed 
only in emission, only in objects with a 
C/O-ratio close to or greater than unity 



and only in cool dusty environments, i.e. 
environments with dust temperatures be- 
low 250 K. Among them are the already 
mentioned C-rich PPNe (about half of 
them show the feature), two Planetary 
Nebulae (Hony, Waters & Tielens 2001) 
with central stars of the Wolf-Rayet type 
and, according to Volk, Xiong & Kwok 
(2000), two extreme carbon stars (IRAS 
21318+5631, IRAS 06582+1507). Clement 
et al. (2004), however, re-analyzing the 
ISO-spectrum of IRAS 21318+5631, find 
no evidence for a feature at 20 /im, and in 
the case of IRAS 06582+1507, Volk, Xiong 
& Kwok (2000) also remark that its reality 
is doubtful. 

It should be noted that an emission 
band with a position and width quite com- 
parable to the '21' /im band is detected in 
the spectra of stars on the early Asymp- 
totic Giant Branch (AGB), i.e. at the 
very opposite end of AGB evolution: the 
19.5 /im feature (Posch et al. 2002). Its 
peak position amounts to 19.4-19.6 /im, its 
profile is symmetric and its full width at 
half maximum (FWHM) is close to 3/im. 
The dust shells of the mostly semiregu- 
lar variable stars in which this feature is 
observed are generally assumed to consist 
of refractory oxides, but not at all of car- 
bonaceous dust grains. A sohd solution of 
FeO and MgO has been proposed as fea- 
ture carrier in this case (see also Ferrarotti 
& Gail (2003)). 

Coming back to the PPNe exhibiting 
the 21 /im band in its 'canonical' form, let 
us consider the spectrum and the derived 
properties of its best-known source, IRAS 
07134+1005 (=SAO 96709 or HD 56126, 
see Fig. 1). The F5 I central star of this 
PPN has a surface temperature of about 
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Fig. 1.— ISO-SWS spectrum of the PPN 
IRAS 07134+1005 at a resolution A/A A 
= 600. At the bottom, the profile of the 
'21'yum band is shown, both as derived 
from the ISO-SWS spectrum HD 56126 
(sohd fine) and from an ISO-SWS spec- 
trum of SAO 34504 (open triangles) by 
subtraction of a polynomial (dotted line 
overplotted to the ISO spectrum at the 
top). 

7000 K according to Hrivnak, Volk & Kwok 
(2000); other PPNe sources of the same 
feature have comparable effective temper- 
atures (5500 to 8000 K). The circumstel- 
lar envelopes of these objects are all op- 
tically thin at infrared wavelengths. As 
for unusual properties of the central stars, 
van Winckel & Reyniers (2000) found them 
to be enriched in s-process elements and 
metal-poor. Apart from the peak due to 
the radiation of the central star at visual 
wavelengths, several infrared bands arising 
from dust grains are present in the spec- 
trum of HD 56126, namely at 8, 11.3, 12.1, 
and 30 ^m. They arc partly unidentified, 
partly attributed to SiC (11.3-broad com- 
ponent), PAHs (11.3-narrow component, 
12+) and magnesium sulfide (30 /um), re- 
spectively. No sihcate features are seen in 
any of the sources of the '21' //m band. 



We derived the profile of the '21' //m 
band using the polynomial subtraction 
method which has been applied by Posch 
et al. (1999) for the derivation of the 13 iim 
band profile. The residual emission of HD 
56126, based on an ISO-SWS-spectrum 
reduced with the ISO Spectral Analysis 
Package (ISAP), is shown at the bottom 
of Fig. 1. For comparison, the scaled resid- 
ual emission of SAO 34504 - calculated by 
the same method - is shown in the bot- 
tom part of Fig. 1 as well (open triangles). 
The close coincidence of the two emission 
profiles is noteworthy; together with the 
result by Volk, Kwok & Hrivnak (1999) 
that the feature profile is pretty much the 
same in all sources, this implies that the 
spectrum of HD 56126 can be considered 
representative and will be so considered 
in the following. The most conspicuous 
properties of the feature profile are a steep 
rise in its blue wing, a maximum position 
of 20.1 ^m and a long 'tail' extending to 
23-24 ^m. 

In contrast to stars on the early AGB, 
bright PPNe have spatial extensions which 
are large enough to resolve them even at 
infrared wavelengths. Kwok, Volk & Hriv- 
nak (2002) and Miyata et al. (2003) tried to 
answer the question whether the different 
broad infrared features of PPNe have their 
origin in different or coincident regions on 
the basis of imaging through broad fil- 
ters. They came to the conclusion that 
the 11.3 um SiC, the 12+ /im plateau emis- 
sion, and the '21'yum emission originate 
(approximately) in the same region, which 
seems to indicate that similar physical and 
chemical conditions are required for the 
formation and/or survival of the respective 
feature carrier. 
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3. Silicon CEirbide, silicon dioxide 
and iron monoxide 

3.1. Silicon ceirbide with impurities 
or coatings 

As mentioned in §1, Speck & Hofmeister 
(2004) suggested silicon carbide with im- 
purities as carrier both of the 11.3 and the 
'21'yum spectral features. One of their ar- 
guments in favour of this hypothesis is the 
spatial co-location of both emission bands 
in some of the objects exhibiting them. 
Obviously, this argument is applicable not 
only to N-doped and C-doped SiC, but also 
to any other particular SiC grain popula- 
tion like e.g. oxidized SiC grains (see be- 
low). 

In the following, when using the term 
'silicon carbide' or 'SiC, we will always re- 
fer to the cubic (beta) modification of this 
polymorphous solid. Due to its symme- 
try, /3-SiC has only one predicted infrared 
active mode, which is located at 796 cm~^ 
(position of the maximum of e"), produc- 
ing, in the case of small spherical particles, 
a maximum of its absorption efficiency fac- 
tor at 935 cm~^ or 10.7 /um. (For elon- 
gated ellipsoidal particles, this maximum 
of the absorption efficiency is shifted to- 
wards longer wavelengths: see Mutschke 
et al. 1999 and below.) As already ex- 
plained, secondary emission bands of SiC 
can occur in the case of a non-ideal crys- 
tal structure, e.g. due to carbon or nitro- 
gen inclusions. The case for N-doped SiC, 
for which Speck & Hofmeister refer to Sut- 
trop et al. (1992), is problematic, how- 
ever, since Suttrop et al. did not only de- 
tect a 20.9 ^m (478 cm~^) band, but also 
stronger additional bands in N-doped SiC. 
These additional bands (e.g. at 667 cm~^ 



(15 //m) and 580cm~^ (17.2//m)) are not 

seen in the sources of the astronomical 
'21' /im feature. Moreover, the whole se- 
ries of nitrogen-induced secondary features 
of SiC is weaker than the SiC main feature 
at least by a factor of 1000. Therefore, 
very low temperatures (below 80 K) would 
be necessary to reverse this ratio of band 
strenghts. According to Hrivnak, Volk & 
Kwok (2000), the dust temperatures at the 
inner radii of the dust shells of the '21' iim 
sources are between 165 and 220 K; due 
to the relatively small ratio between outer 
and inner shell radius, the average tem- 
perature of the dust grains producing the 
'21' //m will hardly be below 80 K even for 
very transparent dust species (see below. 
Fig. 11). 

As for C-doped SiC, there is presently 
too little information at hand about its po- 
tential secondary bands. As long as no 
studies on analytically well-characterized 
SiC samples with carbon impurities are 
available, it remains unclear whether or 
not the properties of N-doped SiC can be 
extrapolated to C-doped SiC. 

Interestingly, there is yet another pos- 
sibility to produce an 11 plus a '21' /xm 
band by, so to speak, 'SiC plus something': 
namely by coating SiC grains with SiO^:. 

Recently performed experiments (Cle- 
ment et al. 2003) demonstrate that pure 
SiC nanoparticles, under the influence of 
the laboratory atmosphere, are quite eas- 
ily oxidized at their surfaces. This ef- 
fect is shown in Fig. 2, where we com- 
pare the transmittance of Si02 grains with 
the transmittance of partially oxidized SiC 
nanoparticles (both embedded in KBr). 
The method for the production of SiC 
nanoparticles by laser pyrolysis and the in- 
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Fig. 2. — Comparison of the transmis- 
sion spectrum of partially oxidized nano- 
SiC particles with that of amorphous Si02 
grains. 
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Fig. 4. — Commercial /9-SiC (dotted line) 
was heated in air at 1050° C (dashed) 
and has afterwards been treated with HF, 
thereby returning from an almost com- 
pletely oxidized state (carrying a 9 and 
a strong 21-22 /xm absorption band) to its 
original state characterized by the well- 
known llyum band. The additional broad 
7 fj,m and 14 /im absorptions in SiC are due 
to plasmon-phonon coupling, see Mutschke 
et al. (1999) for details. 



Fig. 3. — Core-mantle particle consisting 
of SiC and amorphous silicon oxide. The 
silicon carbide core is identified by its char- 
acteristic lattice fringes. 
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Fig. 5. — Volume normalized absorption 
efficiency Qabs/a of core- mantle particles 
composed of an SiC core and an amor- 
phous Si02 mantle, the size of the whole 
(spherical) grains being small compared to 
the wavelength. 



frared spectra of these grains have been 
discussed extensively by Clement et al. 
(2003). They also report the detection 
of an amorphous silicon oxide mantle by 
means of transmission electron microscopy 
(see Fig. 3). The large surface-to- volume 
ratio of the SiC core favours the formation 
and growth of a silicon oxide mantle, which 
can be a few nanometers thick and thus 
reach a considerable volume fraction of the 
particulate. Further evidence for the for- 
mation of Si02 on SiC surfaces, has been 
delivered - even though for significantly 
higher temperatures - by Mendybaev et al. 
(2002) and Tang & Bando (2003). 

Depending on the particle size and 
shape as well as agglomeration effects, 
SiC nanoparticles show one strong IR res- 
onance band, located at 10.8-12.4 yum. 
However, if partially oxidized at their sur- 
faces, the SiC nanograins produce addi- 
tional infrared bands at 9 and 21-22 //m. 
These bands can be assigned to stretching 
and bending vibrations of the Si04 tetra- 
hedra. The minor 12.5 (im band of Si02 is 
not detected in the spectrum of partially 
oxidized SiC nanoparticles, since it almost 
coincides with the SiC main band shifted 
to longer wavelengths by shape effects. 

An additional proof for the oxidation of 
SiC is delivered by treating the oxidized 
grains with HF, which quickly reduces the 
oxidized grains, leaving only the pure sili- 
con carbide cores. The spectroscopic result 
of this process is depicted in Fig. 4. 

Not only experimentally, but also by 
calculating the absorption efficiency of 
core-mantle-particles composed of SiC and 
amorphous Si02, absorption spectra with 
the principal structure of those shown in 
Figs. 2 and 4 can be derived. For this pur- 
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pose, the optical constants of amorphous 
Si02 compiled in the first volume of Palik 
(1985-98) have been used. We performed 
a Lorentz fit to these optical constants in 
order to analyse the underlying oscillator 
parameters. The damping constants of 
the SiC oscillator has been set to 50cm~^ 
and that of the 456 cm~^ oscillator of Si02 
to 42 (instead of 27), assuming that the 
nanocrystalline structure is comparable to 
that of defect-bearing polycrystaUine films 
(see Mutschke et al. (1999), Sect. 2.3). As 
for the 456 cm~^ Si02 oscillator, a larger 
damping constant than that correspond- 
ing to the Pahk data is justifyable also 
on the basis of the spectra published by 
Lehmann, Schumann & Hiibner (1984). 
These authors studied the transmittance 
of thin SiOx films (with x<2) and found 
a strong broadening of the 20-21 //m band 
with decreasing x. This result is of po- 
tential relevance for the coating of SiC by 
silicon oxide in carbon-rich circumstellar 
shells, since in these environments, the for- 
mation of oxygen-deficient Si02 is very well 
conceivable. For a summary of all Lorentz 
parameters used, see Tab. 1. 

Fig. 5 shows the absorption efficiency 
factor Qabs/a of spherical core- mantle- 
particles composed of SiC and amorphous 
Si02 with core volume fractions ranging 
from 40% to 60%. Three prominent max- 
ima of Qabs are discernible: a peak at 
8.3 /im, which originates from the Si02 
mantle; the sharp 11 /um SiC resonance; 
and, finally, a 20 //m peak, with a main 
peak and a shoulder at the long- wavelength 
side, originating again from the Si02 man- 
tle. Peak splitting is typical for bands of 
mantle materials. It arises from the exis- 
tence of two interfaces with core and am- 



Table 1: Frequencies (TOj), osciUator 
strengths {flj) and damping constants (7j) 
of SiC and Si02. We adopted the values 
Coc = 6.49 for SiC and = 2.2 for amor- 
phous Si02. 
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bient medium. All the bands are sharper 
than those of the transmission spectrum 
of partially oxidized SiC shown in Fig. 2. 
This difference in the bandwidth is primar- 
ily due to shape effects. 

Looking at Fig. 5, one may assume that 
the depicted pattern of Cabs/V maxima is 
incompatible with the spectra of astronom- 
ical objects hke HD 56126, since the 20 /im 
band, in this figure, is only a minor peak 
compared to those at 8.2 and 11 /im. How- 
ever, such a conclusion would be premature 
given that the dust in the sources of the 
'21' //m feature is rather cold, as has been 
highlighted in Sect. 2. For dust particles 
with temperatures below 150 K, all fea- 
tures at wavelenghts smaller than ~ 20 /im 
are strongly attenuated, because they tend 
to be located in the Wien regime of the 
corresponding Planck functions. This ef- 
fect, which was already noted by Speck & 
Hofmeister (2004), is illustrated in Fig. 6. 

It is evident from this figure that cold 
core-mantle-particles composed of SiC and 
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Fig. 6. — Products of the absorption ef- 
ficiency for SiC-Si02 corc-mantlc particles 
(core volume fraction = 0.5) with two dif- 
ferent Planck functions, corresponding to 
dust temperatures of 160 K and 100 K. In 
grey colour, the case of eUipsoidal parti- 
cles (shape distribution according to Os- 
senkopf, Henning & Mutschke 1992) is 
overplotted. 
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Fig. 7. — Normalized emission spectrum 
of a dust shell composed of SiC-Si02 core- 
mantle grains. The shell is supposed to fol- 
low a 1/r^ density distribution and a tem- 
perature distribution resulting from the 
condition of radiative equilibrium (see ap- 
pendix) and a certain geometry of the shell 
of HD 56126 as infered from observational 
data (see below, sect. 3.3). The residual 
emission of HD 56126 is shown for com- 
parison. 
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amorphous Si02 do produce a strong 

20 /im feature. There is a rather good 
agreement between the observed and the 
calculated band positions (see also Tab. 
3 for a quantitative assessment), and the 
weakness of the observed 11 /im feature can 
be explained by dust temperatures signif- 
icantly lower than 160 K on the average. 
There remains a problem, however, with 
the width and shape of the 20 /im band 
produced by Si02. For spherical particles, 
the FWHM of the 20 /im SiOa emissivity 
maximum amounts to 1.9 /im, compared 
to an observed FWHM of 2.2-2.3 /xm. Fur- 
thermore, the shoulder at the 'red' wing of 
the 20yum Si02 feature has not been seen 
in the astronomical '21' /xm band. The 
bandwidth discrepancy can be overcome 
by a CDE; the shoulder at the 'red' wing 
of the band, however, does not vanish in 
the CDE case (see Fig. 7); nor does it van- 
ish by any particular (realistic) choice of 
the dust temperature or dust temperature 
distribution. 

3.2. Oxidized silicon grains 

Instead of growing on top of silicon car- 
bide grains, Si02 may equally well form 
mantles on top of pure silicon particles. 
Li & Draine (2002) were the first to exa- 
mine the spectra emerging from silicon 
grains with oxide coatings, focussing on 
the case of stochastically heated nanoparti- 
cles. Smith & Witt (2002) proposed silicon 
nanoparticles as carriers of the so-called 
extended red emission observed in circum- 
stellar environments; however, this identi- 
fication is still a matter of debate. 

The principal structure of the IR spec- 
trum of Si-Si02 core-mantle-grains is simi- 
lar to that of SiC-Si02 core-mantle-grains. 



except for the absence of the 11 fxm band. 

In addition to the 8.2 /im and the 20 fim 
peaks, a rather strong 12.5 /im band ap- 
pears in the Si-Si02 absorption efficiency 
spectrum. The exact position A20 of the 
20 /im resonance depends on the volume 
fraction / of the Si core. A20 increases from 
19.7 to 19.8 nm as / decreases from 0.6 to 
0.4. 



dotted: dust temperature = 160K 
solid: dust temperature = 1 0OK 




X l\im] 



Fig. 8. — Products of the volume normal- 
ized absorption efficiency Cabs/V shown 
in the previous figure with two different 
Planck functions, corresponding to dust 
temperatures of 160 K and 100 K. 

Following the procedure that has been 
explained in the previous subsection, we 
folded the absorption efficiency of a Si- 
Si02 core-mantle-grain with / = 0.5 with 
two different Planck curves B^{Td}. For 
the dust temperature, the two values 
Td = 100 K and Td = 160 K have been cho- 
sen. As for the 20 //m band, the result - 
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plotted in Fig. 8 - strongly resembles the 
case of SiC cores, which is not surprising, 
since the optical constants of the core have 
only a limited influence on the structure 
of those bands which are due primarily to 
the mantle. 

It must be stressed, however, that core- 
mantle particles composed of Si and Si02 
are more unlikely to be the carrier of the 
'21' /xm band than those composed of SiC 
and Si02, since it is questionable whether 
silicon grains can be formed in carbon-rich 
PPNe at all. 

3.3. Cold FeO (wustite) 

As mentioned in §1, iron oxides have 
been considered as carriers of bands in 
the 20-21 yum region by various authors al- 
ready. It seems clear, though, that none 
of the 'higher' iron oxides like Fe203 or 
Fe304 can survive in a reducing environ- 
ment. Iron monoxide (FeO, also called 
wustite^), by contrast, may survive in a 
carbon-rich circumstellar shell. 

According to Duley (1980), the oxida- 
tion of metallic iron is a very efficient pro- 
cess provided that enough O2 molecules are 
present. The degree to which iron parti- 
cles are oxidized depends on the physical 
conditions, inter alia on the temperature. 
The higher iron oxides FC3O4 or FC2O3 be- 
come stable at higher temperatures (see, 
e.g., Gail & Sedlmayr 1998). By contrast, 
as Roberts (1961) and Fehlner & Mott 
(1970) pointed out, several monolayers of 
iron monoxide form very quickly on top 
of Fe surfaces, even at low temperatures, 

^Strictly speaking, wustite is the terrestrial form of 
iron monoxide and has the chemical composition 
Fei_xO with x = 0.05 ...0.16. 



since this process requires a negligible ac- 
tivation energy. In order to end up with 
dust grains composed essentially of FeO 
(instead of a large metallic iron core and 
a tiny FeO mantle) , the (low-temperature- 
) oxidation of very small iron grains is re- 
quired - i.e. grains composed of 10'^ atoms 
at maximum. Only such small iron grains 
- their diameter amounts to ~2 nm - could 
be fully converted into FeO by the process 
mentioned above. The oxidation of iron 
grains consisting of less than 10^ atoms 
is to some extent an artificial assumption; 
less artificial, though, than the constraint 
that the TiC clusters - which have pre- 
viously been considered as carriers of the 
'21' /xm band - should consist of no more 
than 125 atoms to efficiently produce the 
TiC cluster signature. 

Why would FeO, according to our sug- 
gestion, not exist in the vast majority of 
the PNe, at least not with high number 
densities? A plausible explanation is that 
in PNe with their large fractional ioniza- 
tion and high UV photon densities, iron ox- 
ides will be reduced to metaUic iron. While 
the non-detection of the '21' /xm band in 
PNe spectra can thus be accounted for, 
its non-detection in the spectra of most 
of the oxygen-rich as well as carbon-rich 
AGB stars is more difficult to understand 
within the iron grain oxidation scenario 
for PPNe. If in PPNe like HD 56126, 
FeO forms by oxidation of iron, why does 
the same process not take place in C- 
stars, S-stars and in the vast majority of 
M-AGB stars? This may be due to the 
fact that the dust around these stars is, 
on the average, much hotter than dust in 
PPNe and that the sticking probability of 
oxygen on iron decreases as the temper- 



12 



ature increases. The Mg-Fe-oxide forma- 
tion around some M-stars, postulated by 
Posch et al. (2002), must involve an en- 
tirely different mechanism than the low- 
temperature oxidation of iron in PPNe. In- 
deed, Ferrarotti & Gail (2003) predict Mg- 
Fe-oxides via the condensation of the very 
refractory MgO in M-stars. It is notewor- 
thy that HD 56126 and other sources of 
the '21' //m band studied by van Winckel 
& Reyniers (2000) are depleted in metals 
including iron. Other PPNe do not show 
the feature, even though more iron is avail- 
able for FeO formation. There could be 
two reasons for that. One is that the iron 
particles should be small as stated above. 
The second is that FeO formation is essen- 
tially governed by oxygen availability and 
that the iron abundance is therefore only 
crucial if there is enough free oxygen. 

Whereas the knowledge on the Fe- 
Mg-oxide formation in circumstellar en- 
vironments is still scarce, quite detailed 
spectroscopic information on Mg-Fe-oxides 
is available. The optical constants of 
FexMgi_xO samples (x amounting to 0.4 
. . . 1.0) with magnesiowustite structure 
have been measured at room temperature 
by Henning et al. (1995). Subsequently, 
Henning & Mutschke (1997) published low 
temperature data for various cosmic dust 
analogs, i.e. their indices of refraction n(A) 
and absorption k(A) for T = 10, 100, 200 
and 300 K. A major change of the infrared 
spectrum of FexMgi_xO is observed as the 
samples are cooled: namely a sharpening 
of the vibration band located at 18-20 //m 
(depending on x). Fig. 9 shows this effect 
for x = 1.0, i.e. for pure FeO. (The case 
T = 10 K has been omitted, since it does 
sot significantly differ from T = lOOK.) 




0,0- 
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Fig. 9. — Absorption efficiency factors Cabs 
of spherical FeO grains, calculated from 
their optical constants measured at room 
temperature (T = 300 K) and at low tem- 
peratures (T = 200K, T = 100K). The re- 
spective Cabs functions have been multi- 
plied with Planck functions B^, (T) (nor- 
malized at 20 /im) for the appropriate tem- 
peratures. 
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The FWHM of the 20 /im FeO band de- 
creases from 3.6 fim at room temperature 
to 2.4 //m at T = 100 K, while the band po- 
sition is shifted from Amax= 19.9 /xm to 
Amax = 20.1/xm. The latter effect is exclu- 
sively due to the multiplication with the 
Planck functions B,^ (T), which have max- 
ima positions of 17, 25.5 and 51 /xm for 
T = 300, 200 and 100 K, respectively. The 
rather large width of the 20 //m band pro- 
duced by FeO, in combination with the fact 
that for T = 200 K and T = 100 K, the peak 
position is located at the 'blue' wing of 
the corresponding Planck curve, leads to 
an asymmetric profile of the CabsXBi,(T) 
curves - the same kind of asymmetry which 
is found for the astronomical '21' /im band. 

Iron monoxide has a rather high absorp- 
tion efficiency in the visual range of the 
spectrum. Therefore, it is comparatively 
strongly heated in the radiation field of an 
F or G star. An equilibrium temperature 
of less than 100 K is reached only at dis- 
tances from the central star larger than 
2 X 10^ stellar radii. This can be shown 
by calculating the energy balance for the 
case of an equilibrium between the radia- 
tion field of the central star and the dust 
grains in its circumstellar according to the 
scheme summarized in the appendix. 

According to Kwok, Volk & Hrivnak 
(2002), the inner radius Rin of the dust 
shell of IRAS 07134+1005 is located at 
1.2" from the central star, corresponding 
to 1.2 (1.8) X 10~^ pc for an assumed dis- 
tance of 2 (3) kpc (Kwok, Volk & Hrivnak 
(2002) assume a distance of 2 kpc, Meixner 
et al. (2003) assume 3 kpc). In terms of 
the stellar radius corresponding to a lumi- 
nosity of 1400 L0 kpc"^ (Hrivnak, Volk & 
Kwok 2000) and a temperture of 7000 K, 



this amounts to 9760 or 6490 stellar radii 
for d=2 (3) kpc, respectively. 



1000 




10"^ 10' 10* 10' 
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Fig. 10. — Radiative equilibrium tem- 
perature of FeO in the radiation field 
of a star with an effective temperature 
of 7000 K (solid line). For comparison, 
the radial temperature profile dust with 
wavelength-independent optical constants 
is also shown (dashed line; sec also the ap- 
pendix to this paper). Short-dotted hues 
dehmit the r-range of the brightest part of 
IRAS 07134+1005. 

In Fig. 10, the position of Rin/Reg = 6490 
has been highhghted by a dotted verti- 
cal line. The equilibrium temperature of 
FeO corresponding to this distance from 
the central star amounts to 155 K. This 
is in good agreement with the results of 
the model calculations by Hrivnak, Volk & 
Kwok (2000), who derive an inner temper- 
ature of the dust shell of 165 K by radiative 
transfer calculations aimed at a fit of the 
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spectral energy distribution. 
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Fig. 11. — Equilibrium temperatures of 
amorphous carbon, FeO and core-mantle 
particles composed of SiC and Si02 in the 
radiation field of an F5 star. The do- 
main from which presumably most of the 
IR emission of IRAS 07134+1005 comes is 
indicated by dotted lines. 

In terms of the ratio between outer 
and inner shell radius, the dust shells of 
the PPNe showing the '21' //m feature are 
rather thin. For IRAS 07134+1005, the 
observations of Kwok, Volk & Hrivnak 
(2002) and Miyata et al. (2003) yield an 
outer shell radius not much larger than 
twice the inner shell radius. This im- 
plies that the dust temperature range, 
for a r~'''^-like decrease of T(r),'^ is rather 
small. For FeO, we derive a dust tem- 

^See the appendix for a comment on the meaning 
of the T(r) oc r~° '' decrease. 



perature at 2 Rin of 120 K, meaning that 
120K<Td<155K (again for d=3kpc). For 
other dust species which are generally as- 
sumed to be present in carbon-rich PPNe, 
narrow T^ ranges are found as well. Amor- 
phous carbon is warmer than FeO due to 
the monotonous increase of its absorp- 
tion efficiency factor from the IR to the 
UV: we find 153K<Td<200K under the 
same conditions as above; core-mantle par- 
ticles composed of SiC and Si02, unless 
heavily N-doped, remain cooler than FeO 
(104K<Td<126K;^ sec Fig. 11). 

These estimates make it possible to cal- 
culate the emission spectrum of an opti- 
cally thin dust shell with a homogeneous 
density distribution and a temperature 
range of 120 K<Td<155 K. If we are only 
interested in the profile of the emerging 
spectral energy distribution, it suffices to 
calculate a mean of blackbody functions 
B^(T) for 120K<Td<155K and to mul- 
tiply the result with Cabs(Tm)- The re- 
sulting normalized spectrum is shown in 
Fig. 12. Both with respect to feature posi- 
tion and width, the emission spectrum of 
FeO calculated for the temperature range 
indicated above reproduces the observed 
'21' fim band profile fairly accurately, even 
though no ad hoc assumptions concerning 
the dust temperature have been made. 

For a CDE, quite different spectra of 
FeO emerge. For a CDE according to Os- 
senkopf, Henning & Mathis (1992) ('CDE- 
Osk'; it is characterized by a maximum 
probability for particle shapes close to the 
spherical shape), the position of the FeO 
absorption efficiency maximum is shifted 

^These values refer to a core volume fraction of 0.5; 
the larger the (rather transparent) Si02 mantle, 
the lower the dust temperature will be. 
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Fig. 12. — Normalized emission spectrum 
of an FeO dust shell (solid: spherical parti- 
cles, dotted: CDE according to Ossenkopf 
et al. (1992)) with a temperature distribu- 
tion corresponding to the range highlighted 
in Fig. 11 and a 1/r^ density distribution. 
For comparison, the residual emission of 
IRAS 07134+1005 is shown. 



to 20.7-21 iim, and the bandwidth changes 

even more dramatically, i.e. it increases to 
6.5-7.6 /im (see Fig. 12). For a CDE with 
equal probability of all particle shapes, the 
changes in feature position and shape are 
even more dramatic, but this kind of shape 
distribution is unlikely to represent the 
shapes of circumstellar grains. The as- 
sumption that nearly spherically symmet- 
ric grains dominate in circumstellar dust 
shells is more reasonable, hence the 'truth' 
should lie between the 'CDE-Osk' case and 
the case of spherical particles. 

Fig. 12 may give the impression that 
FeO dust can produce a feature compara- 
ble to the astronomical '21' //m emission 
band only under very particular phys- 
ical conditions, e.g. only for T<155K, 
while wc know from Hrivnak, Volk & 
Kwok (2000) that the dust temperatures 
in the sources of the '21' //m band can 
reach 220 K. However, as long as T<250 K, 
Bi^(T) will always have its maximum at a 
wavelength larger than 20 ^m, thus shift- 
ing the 19.9 //m emissivity maximum of 
FeO to a feature position of 20.0-20.1 ^m. 
And indeed, T<250K seems to be fulfilled 
for all known sources of the '21' /im fea- 
ture. In a wider sense, however, it is true 
that FeO can produce a strong 20.1 /xm 
feature only under particular physical con- 
ditions. If the mean dust temperature is 
much smaller than 100 K, too little en- 
ergy will be contained in the 20 /im FeO 
band to produce a feature as strong as 
that observed in some PPNe. If, by con- 
trast, T>250K, the feature will become 
too broad and to much shifted to the 'blue' 
as to account for the astronomical '21' /im 
band.^ 

^Henning & Mutschke (1997) point out that only 
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3.4. Remcirk on the elemental abun- 
dances required to account for 
the '21' /xm band 

Li (2003) made a very useful estimate 
of tlie mass of dust whicli is required to 
account for the radiative energy contained 
in the '21' /im band in the case of IRAS 
07134+1005. He pointed out that the 
Kramers-Kronig equation relates the in- 
tegral over the absorption efficiency of a 
grain to its volume and the orientational 
averaged polarizability F. Since for all ma- 
terials with a static dielectric constant >4 
and for moderately elongated grain shapes, 
the latter quantity has a value between 0.5 
and 1 (Purcell 1969, Fig. 1), there is an 
upper limit to the integrated absorption 
efficiency per unit volume, which is more 
or less the same for the dust species con- 
sidered here. Consequently, there is also 
a general upper limit to the energy emit- 
ted in the 21 //m band per unit particle 
mass or, vice versa, a minimum total grain 
volume for a certain total energy emitted 
through this band. Hence, the minimum 
mass derived by Li can be adopted to other 
materials just by scaling with the mass 
density of the grain material. This leads 
to the following conclusion concerning the 
abundance problem: While there is, ac- 
cording to Li (2003), about 100 times too 
little amount of titanium to account for the 
'21' //m band by assigning it to TiC, there 
is 2.5 times more Fe than would be needed 
to account for the '21' micron band with 
FeO (Fe is 250 times more abundant than 
Ti in HD 56126). Sulfur is 312 times more 

a limited number of cosmic dust analogs have op- 
tical constants which change as strongly at low 
temperatures as those of FeO (and FeS). 



abundant than Ti in the same star, silicon 
1225 times, carbon 34500 times and oxygen 
36000 times. Hence, compounds of iron, 
sulfur, silicon, carbon and oxygen cannot 
be excluded as carriers of the '21' yum band 
in HD 56126 and objects with comparable 
elemental abundances. 

4. Discarded feature carriers 

For the sake of completeness, two po- 
tential carriers of the '21' /im band which 
have been considered as promising candi- 
dates (but should no longer be considered) 
will be discussed: silicon disulfide and ti- 
tanium carbide. In the case of TiC, insur- 
mountable problems with the abundance of 
titanium exist, whereas SiS2 can be ruled 
out by spectroscopic arguments. 

4.1. Silicon disulfide 

As mentioned in §1, the formation of 
sulfides is very likely in carbon-rich cir- 
cumstellar environments. Indeed, the 
26 iJ,m feature, observed in many sources of 
the '21' //m feature, is currently ascribed 
to magnesium sulfide (MgS, Hony et al. 
2003). 

However, from a spectroscopic point 
of view, there are problems with ascrib- 
ing the 20.1 //m band to SiS2. The main 
point is that this material has two in- 
frared bands, of which the second, lo- 
cated at 16.8 /im for spherical particles 
(see Fig. 13), has no counterpart in the 
astronomical spectra. Furthermore, the 
position and FWHM of the main SiS2 
band, which amout to 19.8 and 1.6 /im, re- 
spectively, do not agree with the observed 
characteristics of the 20.1 /xm band. For 
a CDE, this discrepancy becomes smaller 
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(Amax = 20.2 /xm, FWHM=1.8/xm), but 
the problem with the unobserved sec- 
ondary feature at 16.8 /im becomes even 
worse, given that the ratio between ~17 
and ~20//m band increases for a CDE. 
Very low dust temperatures (significantly 
below 100 K) would be required to make 
the ~17/im band negligible in strength 
compared to the ~20 /im SiS2 band. 

The visual and near-infrared optical 
constants of SiS2 have not yet been de- 
rived, such that it is not possible to cal- 
culate T(r) for SiS2 based on its radiative 
energy balance. However, it is hardly con- 
ceivable that SiS2 is so transparent in the 
visual range as to remain much colder than 
100 K over the whole range indicated in 
Fig. 11. 

4.2. Titanium carbide 

Even though recent papers by Li (2003), 
Chigai et al. (2003) exclude TiC nanopar- 
ticles as the carrier of the '21' //m band, 
we append a discussion of TiC's infrared 
properties here, focussing on the relation 
between bulk TiC and TiC clusters. 

Macroscopic TiC has metallic proper- 
ties, with an absorption index k which is 
larger than the refractive index n through- 
out the whole infrared range. For A>13 /im, 
both n and k are larger than 10. At about 
19 /xm, a resonance peak occurs, which is 
rather weak, however, with a relative in- 
crease of n, k and Cabs by less than 10%. 

Small TiC clusters, e.g. those composed 
of 4x4x4 atoms, behave differently from 
bulk matter insofar as the mobility of the 
cluster electrons is reduced. The lattice 
resonance is then no longer a small ad- 
ditional term in the dielectric function, 
but the dominant contribution at infrared 




Fig. 13. — Normalized emission spectrum 
of spherical SiS2 particles for two differ- 
ent dust temperatures, 100 and 160 K. The 
16.8 nm secondary emission band, which is 
not observed in the sources of the '21' //m 
band. 



18 



frequencies. Therefore, a comparatively 
strong 20 /xm emission band arises from the 
clusters; its exact position, however, de- 
pends on the cluster size. Chigai et al. 
(2003) derived the following Lorentz oscil- 
lator parameters for the small TiC clusters 
produced by von Helden et al. (2000): a 
resonance frequency of 488 cm^^, a damp- 
ing constant of 50.2 cm~^ and an oscillator 
strength of 187 cm~^. For spherical parti- 
cles, this translates into a maximum posi- 
tion of the absorption efficiency factor of 
20.0 ^m and an FWHM of 2.0 fim. For a 
CDE, no shift of the band position is found 
in this case (due to the non-negative values 
of the dielectric function of TiC according 
to Chigai et al. (2003)), but only an in- 
crease of the FWHM to 2.2 /im (for a CDE 
according to Ossenkopf, Henning & Mathis 
1992). 

Fig. 14 shows the normahzed emission of 
a small TiC cluster (at Td=160K), the nor- 
malized dust emission arising from macro- 
scopic dust grains (at two different tem- 
peratures Td =160K and T(j=120K) and, 
again for comparison, the profile of the as- 
tronomical '21' yum feature as derived from 
the ISO-spectrum of IRAS 07134+1005. 
Since the maximum of the absorption ef- 
ficiency factor of TiC is not evident from 
Fig. 14 due to the normalization, it should 
be added that Cabs/V (TiC) does not ex- 
ceed 0.6 /im~^. (Depending on the value 
of eoo, Cabs/V may also be substantially 
smaller). By contrast, the ~20 /im peak 
absorption efficiency of cold (100 K) FeO 
as well as that of amorphous Si02 amounts 
to 1.6 jim^^. 
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Fig. 14. — Normalized emission spectrum 
of bulk TiC at temperatures of 120 K 
(dashed) and 160 K (sohd) and of TiC 
nanoparticles at T = 160 K (dotted) com- 
pared to the residual emission of IRAS 
07134+1005. 



5. Conclusions 

Tables 2 and 3 summarize the spectro- 
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scopic information on the potential carri- 
ers of the '21' /im feature in two different 
respects. Tab. 2 gives an overview of the 
positions of the absorption efficiency max- 
ima and of the corresponding bandwidths. 
Tab. 3, on the other hand, refers to normal- 
ized dust emission spectra like those shown 
in Figs. 7, 12, 13 and 14, hence to products 
of the absorption efficiency with blackbody 
functions. These products are compared 
with the repeatedly shown profile of the 
'21' /im band: we present the sums (hence- 
forth: of the squares of the differ- 
ences between the 'profile' and the 'dust 
emission' values, for wavelengths between 
17 and 25 yum (outside this range, the ISO 
profile is characterized by comparatively 
large noise). The 'temperature grid' has 
been chosen such as to cover the range be- 
tween 160 K and 70 K in steps of 10 K. The 
bold numbers indicate the respective min- 
ima of S(52. For SiC+SiOs-sph, TiC-sph 
and FeO-sph, the minimum of occurs 
at blackbody temperatures between 80 and 
90 K. This is due to the asymmetric shape 
of the profile of the '21' //m band. On the 
basis of a symmetric aborption efficiency 
profile, this asymmetric profile is best re- 
produced by multiplication with a black- 
body of 80-90 K. If, however, the absorp- 
tion efficiency of the feature carrier has an 
intrinsically asymmetric profile (as for el- 
lipsoidal SiC-|-Si02-core-mantle grains), it 
is not necessary to assume such a low dust 
temperature (120 K are sufficient in this 
case). The values for SiS2 are partic- 
ular insofar as they continuously decrease 
with decreasing temperature, down to less 
than 70 K. This has a very simple rea- 
son: The 16.8 /jLUi secondary band of SiS2 
(see Fig. 13), which has no counterpart in 
the astronomical spectra, is the more sup- 



pressed the lower the temperature. As re- 
gards the absolute values of they are 
smallest for TiC clusters and second small- 
est for ellipsoidal SiC-|-Si02 core-mantle- 
grains. 

Despite the good spectroscopic fits ob- 
tainable with TiC, cold SiC-|-Si02 and cold 
FeO, we come to the conclusion, based on 
the properties of the sources of the '21' yum 
band, that none of the hitherto discussed 
candidate carriers is a convincing feature 
carrier candidate in all relevant respects. 
Even though there are oxides which may 
appear as likely band carriers from a spcc- 
troscopical point of view, it is difficult to 
conceive of a scenario accounting for their 
formation in carbon-rich PPNe. On the 
other hand, dust species which are usually 
predicted to form in carbon-rich environ- 
ments, like carbon, carbides and sulfides 
do either not match the observed '21' //m 
band profile (e. g. 8182), or even their max- 
imum abundance is too small for emitting 
the energy contained - at least in some 
cases - in the '21' /xm feature (which rules 
out TiC). As for the individual groups of 
feature carriers, our arguments are the fol- 
lowing ones: 

1) Dust and macromolecular species 
which produce a 20-22 fim emission by 
transitions involving hcteroatoms in car- 
bonaceous networks (e.g. Papoular 2000), 
are not likely to be the feature carrier 
due to the uniqueness of the '21' /xm band 
profile on the one hand and the non- 
uniqueness of the heteroatom-caused emis- 
sion bands on the other hand. 

2) There do exist oxide dust species 
which can produce a 20.1 yum emission 
band very effectively (if the question of 
their formation in C-rich environments is 
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Table 2: Overview of the most promising carriers of the '21'yum feature which have been 
discussed in the present paper. Note: The observed position and FWHM of the '21' //m 
band are 20.1 /xm and 2.3 /um, respectively, 'sph' means spherical particle shape, 'CD&Osk' 
designates a continuous distribution of ellipsoids according to Ossenkopf, Henning & Mathis 
(1992). The A (Cabs,max) values do not directly represent feature positions, since these result 
from the product CabsXBi/ (Tj); for dust temperatures close to 150 K, this leads to a shift 
towards the 'red' of 0.1-0.2 //m. 



candidate carrier particle shape 


A (Cabs.max) [/^l^] 


FWHM (20) [i^m] 


SiC+SiOs 


sph.*^"^ 


8.1-8.3, 11.0-11.1, 12.5, 19.7-19.8 


1.9 




CDE-Osk(«) 


8.2-8.4, 11.1-11.2, 19.8-20.0 


2.1-2.6 


Si+SiOs^") 


sph. 


8.0, 12.5, 19.7-19.8 


1.8 


FeO(^) 


sph. 


19.9 


2.3 




CDE-Osk 


20.7 


6.5 


SiS2 


sph. 


16.8, 19.8 


1.6 




CDE-Osk 


16.8, 20.2 


1.8 


TiC clusters 


sph 


20.0 


2.0 




CDE-Osk 


20.0 


2.2 


TiC grains 


sph 


19.2 


-(c) 



core-mantle- particles with SiC core volume fractions between 0.4 and 0.6 
(''^at an average temperature of ~ 140 K 
('^^broad bump, FWHM not defined 



21 



Tabic 3: Selected potential carriers of the '21' //m feature: quality of fits to the profile of 
the '21' //m band for different mean dust temperatures (between 130 K and 70 K). As 
an indicator for the quality of an individual fit, the sum of the squares of the differences 
'normalized 21 fim fiux' minus 'Cabs x (T^)' has been calculated for 17 fim < A < 25 /im. 
This sum corresponds to the usual 'x^' insofar as the quantities subtracted from each other 
are already normalized. 



[K] 


SiC+SiOa-sph^"^ 


SiC+Si02-CDE^«^ 


FeO-sph 


TiC-sph 


SiS2-sph 


160 


5.98 


3.10 


6.40 


3.54 


18.3 


150 


5.66 


2.92 


6.03 


3.29 


17.8 


140 


5.31 


2.80 


5.62 


3.02 


17.4 


130 


4.94 


2.71 


5.20 


2.72 


16.9 


120 


4.55 


2.67 


4.78 


2.41 


16.4 


110 


4.14 


2.74 


4.37 


2.08 


15.8 


100 


3.76 


3.00 


4.03 


1.76 


15.1 


90 


3.47 


3.63 


3.87 


1.50 


14.3 


80 


3.42 


4.99 


4.13 


1.40 


13.4 


70 


4.02 


7.94 


5.41 


1.74 


12.4 



core-mantle-particles with a SiC core volume fraction of 0.5 



left aside). As wc have shown, this holds 
true both for cold FeO and for Si02. Iron 
monoxide has only one IR active mode, 
centered at 19.9 //m, which shows an ex- 
ceptional narrowing at temperatures in the 
200-100 K range. By contrast, amorphous 
Si02 has two strong IR active modes, one 
of them being centered at 20.5 ^m. Our 
laboratory experiments show that amor- 
phous Si02 easily forms mantles on top 
of SiC. Core-mantle particles composed of 
SiC and Si02, if colder than 150 K, emit 
mainly through the lower energy band of 
Si02. Thus they can in principle account 
for the '21'yum band, even though the re- 
production of the feature profile is not as 
good as in the case of cold FeO. 

But there is a major problem with any 
scenario involving oxides as band carri- 



ers: the difficulty to explain how oxides 
can form or at least survive in carbon-rich 
PPNe. It is hardly conceivable that the 
carrier of the '21'yum band is just a rem- 
nant of the oxygen-rich outflow phase of 
the respective former AGB stars. The in- 
terval between the oxygen-rich phase and 
the carbon-rich PPN phase is too long as 
to account for a co-location, as proven for 
HD 56126, of products of the C-rich and 
the 0-rich phase. The remnants of the O- 
rich phase are expected to be located at 
least one order of magnitude further out 
on a radial scale than the products of the 
C-rich phase. How, then, can either SiC 
be oxidized or FeO form/survive in C-rich 
environments? As for the first process, 
there are several studies showing that com- 
plete transformation of SiC to Si02 occurs 
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at temperatures above 700 K (e.g. Mendy- 
baev et al. 2002, Tang & Bando 2003); pro- 
vided that large defective surfaces as those 
of nanograins are available for the reaction, 
the formation of Si02 mantles on top of 
SiC grains also takes place at room tem- 
perature as our laboratory experiments in- 
dicate (see Fig. 3). But why should this 
Si02 mantle formation take place prefer- 
entially in PPNe with C/0 ratios close to 
or larger than unity? This question, un- 
fortunately, has to remain open here. As 
for the hypothesis that FcO is the carrier 
of the '21' /im band, it would require the 
oxidation of Fe to FeO in cold circumstel- 
lar environments. It seems that metal-poor 
PPNe hke HD 56126 - paradoxically at the 
first glance - provide good conditions for 
FeO formation, since this process requires 
the presence of small metallic iron grains. 
The absence of any feature attributable to 
FeO in the vast majority of the Planetary 
Nebulae, on the other hand, could be ac- 
counted for by the fact that FeO can per- 
sist only in a narrow range of equilibrium 
between oxidation (to higher iron oxides) 
and reduction (to metallic iron). Thus, in 
PNe, FeO is probably reduced to metallic 
iron. 

3) Quantum size effects critically influ- 
ence the optical properties of very small 
grains and clusters. Phonon modes can 
shift when the lattice changes from 'infi- 
nite' to finite size, and conductivity may be 
reduced. For the TiC grainlets proposed 
as carriers, this means that a '21' //m band 
would only occur for clusters and particles 
in a very narrow size range. This makes 
it unlikely that there would be a sufficient 
number of TiC particles in the PPNe to 
account for the observed band strength. 



given that the predicted Ti abundance is 
not large enough and that mantling by 
amorphous carbon will reduce the band 
strength even further. 

4) Any identification or ruling-out of a 
candidate carrier is strengthened if it does 
not only have one, hut two or more ob- 
servable infrared bands. If these bands 
are seperated by >10/im, the secondary 
bands may be suppressed by temperature 
effects. This is likely to occur for the 11 //m 
SiC band expected to arise from the core- 
mantle particles mentioned above. How- 
ever, in the case of SiS2, even tempera- 
tures in the order of 100 K will not be suf- 
ficient to suppress the 16.8 //m secondary 
emission feature. The multiband identi- 
fication problem also sets constraints on 
bands produced by impurities, because the 
impurity-induced bands are usually much 
weaker than the phonon bands of the re- 
spective matrix materials (which provides 
a further argument against doped SiC as 
the carrier of the '21' /im feature). 

We thank Gabriele Born, Jena, for 
preparing our samples. Ingrid Hodous, Vi- 
enna, kindly provided us with the ISO- 
spectra of IRAS 07134+1005 and SAO 
34504, reduced with the ISO spectral anal- 
ysis package ISAP. An anonymous referee 
contributed some substantial arguments 
to the original version of this paper and 
thereby improved it. TP received a DOC 
grant from the Austrian Academy of Sci- 
ences. HM acknowledges support by DFG 
grant Mu 1164/5. 
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A. Appendix: On the temperature profile and spectral energy distribution of 
geometrically thin PPNe dust shells 

In the previous sections, calculations of the radiative equilibrium temperature of the 

sources of the '21' /xm feature have been presented. The temperature profiles in turn have 
been partly used to calculate model spectra. We append a discussion of the method used 
for the dust temperature (Tj) calculation as well as of the influence of on the emerging 
spectral energy distributions (SEDs). 

In the radiative equilibrium case and for optically thin shells (where backwarming effects 
can be neglected), the radial variation of is derived from the following energy balance: 

/•oo 7->2 rco 

/ 7ra^gabs(A)7rSA(A,Teff)— s^dA= / AT:a^Q,,,,{\)TiBx{X,T^) d\ (Al) 

The first term describes the energy absorbed by an individual dust grain with a radius a 
and an absorption efficiency factor (5abs(A) in a stellar radiation field Bx{\ T^^) {T^^ denoting 
the effective stellar temperature); the second term describes the energy emitted by the same 
dust grain. i?cfr is the effective stellar radius, whereas R{T^) is the distance of the grain with 
temperature T^ from the star. While the above integral runs from zero to infinity, numerical 
integration can of course be performed only between a certain lower (Amin) and a certain 
upper limit (Amax)- We chose Amin = 0.05 //m and Amax — 500 //m. It is difficult, however, 
to derive the values of Qahs close to Amin; some of the sets of optical constants have to be 
extrapolated to UV wavelengths since measurements of the UV properties of cosmic dust 
analogues arc still lacking. For FcO, e.g., our measurements extend to 0.2 /zm; beyond this 
wavelength, n and k were assumed to be constant. Any other assumption, and even the 
choice Amin = 0.2 /xm, would not significantly change our results. 

We recall that simple power-law models for Qabs(A) make it possible to solve eq. (Al) 
analytically for T(r). For Qabs(A) oc 1/A - the case corresponding to small dust particles 
with wavelength-independent (practically: weakly wavelength-dependent) optical constants 
n and k -, T(r) decreases with l/r°'^. This solution to eq. (Al) has been plotted in Fig. 10. 
(For the less realistic case Qabs(A) = const., T(r) decreases with l/r°-^.) 

In the main part of our paper, it has been shown that the temperature range corresponding 
to the measured spatial extension of a typical PPN hke HD 56126 is comparatively narrow 
for all the dust species examined. The ratio between the temperature at the outer and the 
temperature at the inner dust shell edge (i.e., its dense observable part) is typically in the 
order of 1.3. This has an important consequence for the calculation of the dust shell's SED. 
Again in the optically thin case, the latter is given by 

/•r=Rout/Rin 

i^d(A) = c X gabs(A) / 5a(A, Td(r)) dr, (A2) 

Jr=l 

where c is a product containing the dust column density, the average geometrical cross 



24 



section of the dust grains, and the surface of the inner dust shell boundary (e.g. Schutte 
& Tielens (1989)). If Kout/^in is indeed in the order of 2 (as for HD 56126 according to 
the observations of Kwok, Volk & Hrivnak (2002) and Miyata et al. (2003)), the result of 
the integral on the right-hand side of eq. (A2) is surprisingly similar to a single black-body 
function for a temperature intermediate between T(i(Rin) and T(j(Rout)- This effect is shown 
in Fig. 15, where we have set c=l and Qabs(-^)=1 for simplicity. 




10 100 

X [|jm] 



Fig. 15. — Spectral energy distributions (Fi,(A)) emerging from a single blackbody (dotted 
line, for T=175K) compared to SEDs emerging from a sum of blackbody contributions as 
the one entering into eq. (A2). For Rout/Rj„ < 2 (sohd and dash-dotted line), the blackbody- 
sum is virtually not distinguishable from a single blackbody with a temperature intermediate 
beween T(Rj„) and T(Rout)- Only for large values of Rout/^in (e.g. values close to 10000, 
as for the dashed line), the blackbody-sum has a profile distinctly different from that of a 
single Planck function. 

It is evident from this figure that the the SED calculated for a geometrically thin cir- 
cumstellar shell (using realistic absorption efficiencies) according to eq. (A2) depends much 
more on the order of magnitude of the dust temperature (i. e. on whether the mean dust 
temperature amounts to 175 or 125 K) than on actually calculating a weighted average of 
Planck functions. 

On the other hand, even a reasonably good estimation of the mean dust temperature 
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Td(r) is difficult without any information on the values of Qabs('^) in the visual and near- 
infrared range. For SiS2 and TiC, e.g., no optical constants for wavelengths smaller than 
10 fim are available. Hence, the SEDs emerging from dust shells composed of these grain 
species can only be calculated by arbitrarily assuming mean T^-values for them. 
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